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Global three— neutrino oscillation analysis of neutrino data 
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Abstract 

A global analysis of the solar, atmospheric and reactor neutrino data is 
presented in terms of three-neutrino oscillations. We include the most re- 
cent solar neutrino rates of Homestake, SAGE, GALLEX and GNO, as well 
as the recent 1117 day Super-Kamiokande data sample, including the recoil 
electron energy spectrum both for day and night periods and we treat in a 
unified way the full parameter space for oscillations, correctly accounting for 
the transition from the matter enhanced (MSW) to the vacuum oscillations 
regime. Likewise, we include in our description conversions with 612 > 7r/4. 
For the atmospheric data we perform our analysis of the contained events and 
the upward- going i^-induced muon fluxes, including the previous data samples 
of Frejus, 1MB, Nusex, and Kamioka experiments as well as the full 71 kton- 
yr (1144 days) Super-Kamiokande data set, the recent 5.1 kton-yr contained 
events of Soudan2 and the results on upgoing muons from the MACRO detec- 
tor. We first present the allowed regions of solar and atmospheric oscillation 
parameters 612, Am^i and ^23; ^"^32! respectively, as a function of ^13. We 
determine the constraints from atmospheric and solar data on the mixing an- 
gle 013, common to solar and atmospheric analyses. The solar limit on ^13, 
although relatively weak, is totally independent on the allowed range of the 
atmospheric mass difference Am'^2- ^^ the other hand the atmospheric data 
analysis indicates an important complementarity with the reactor limits al- 
lowing for a stronger constraint on the allowed value of ^13. We also obtain the 
allowed ranges of parameters from the full five-dimensional combined analysis 
of the solar, atmospheric and reactor data. 
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I. INTRODUCTION 

Super-Kamiokande high statistics data [1-3] indicate that the observed deficit in the /i- 
hke atmospheric events is due to the neutrinos arriving in the detector at large zenith angles, 
strongly suggestive of the z/^ oscillation hypothesis. Similarly, their data on the zenith angle 
dependence and recoil energy spectrum of solar neutrinos [4,5] in combination with the 
results from Homestake [6], SAGE [7], and GALLEX+GNO [8,9] experiments, have put on 
a firm observational basis the long-standing problem of solar neutrinos, strongly indicating 
the need for v^ conversions. 

Altogether, the solar and atmospheric neutrino anomalies [1-5,10-13] constitute the only 
solid present-day evidence for physics beyond the Standard Model (SM). It is clear that 
the minimum joint description of both anomalies requires neutrino conversions amongst all 
three known neutrinos. In the simplest case of oscillations the latter are determined by 
the structure of the lepton mixing matrix [14], which, in addition to the Dirac-type phase 
analogous to that of the quark sector contains two physical [15] phases associated to the 
Majorana character of neutrinos. CP conservation implies that lepton phases are either zero 
or TT [16]. For our following description it will be correct and sufficient to set all three phases 
to zero. In this case the mixing matrix can be conveniently chosen in the form [17] 



' C13C12 S12C13 Si3 \ 



R 



(1) 



^Si2C23 — S23S13C12 C23C12 — S23S13S12 S23C13 
\ S23S12 - S13C23C12 -S23C12 - S13S12C23 C23C13 J 

With this the parameter set relevant for the joint study of solar and atmospheric conversions 
becomes five-dimensional: 
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where all mixing angles are assumed to lie in the full range from [0,7r/2]. 

In this paper we present a global analysis of the data on solar, atmospheric and reactor 
neutrinos, in terms of three-family neutrino oscillations. There are several three-neutrino 



oscillation analyses in the literature which either include solar [18,19] or atmospheric neutrino 
data [20]. Joint studies were also performed, but without including the most recent and 
precise Super-Kamiokande data [21]. This work updates and combines all these results in 
a unique comprehensive analysis. It is known that in the case Aml2 ^ A^^i2) fo^ ^i3 = 
the atmospheric and solar neutrino oscillations decouple in two two-neutrino oscillation 
scenarios. In this respect our results also contain as limiting cases the pure two-neutrino 
oscillation scenarios and update previous analyses on atmospheric neutrinos [22,23] and solar 
neutrinos [24,25] (for an updated analysis of on two-neutrino oscillation of solar neutrino 
data see also Ref. [26]). 

We include in our analysis the most recent solar neutrino rates of Homestake [6], 
SAGE [7], GALLEX and GNO [8,9], as well as the recent 1117 day Super-Kamiokande data 
sample [5], including the recoil electron energy spectrum both for day and night periods. 
Treating in a unified way the full parameter space for oscillations, we carefully account for the 
transition from the matter enhanced to the vacuum oscillations regime. Likewise, we include 
MSW conversions in the second octant (dark-side), characterized by 612 > 7r/4 [18,19,27]. 
As for atmospheric neutrinos we include in our analysis all the contained events as well as 
the upward-going neutrino-induced muon fluxes, including both the previous data samples 
of Frejus [28], 1MB [10], Nusex [29] and Kamioka experiments [11] and the full 71 kton-yr 
Super-Kamiokande data set [3], the recent 5.1 kton-yr contained events of Soudan2 [12] 
and the results on upgoing muons from the MACRO detector [13]. We also determine the 
constraints implied by the CHOOZ reactor experiments [30]. 

From the required hierarchy in the splittings Am^^^ ^ AttIq indicated by the solutions to 
the solar and atmospheric neutrino anomalies (which indeed we justify a posteriori) it follows 
that the analyses of solar data constrain three of the five independent oscillation parameters, 
namely, Am2i,^i2 and ^13. Conversely, the atmospheric data analysis restricts Am|2) ^23 
and ^13, the latter being the only parameter common to both and which may potentially 
allow for some mutual influence. In our global approach we will statistically combine these 
solar and atmospheric limits on ^13. We also compare this bound and combine them with 
the direct limit on ^13 which follows from reactor experiments. While the solar limits on ^13 
are relatively weak, our atmospheric data analysis indicates an important complementarity 
between the atmospheric and the reactor limits. 

The outline of the paper is the following: in Sec. II we review the theoretical calculation 



of the conversion probabilities for solar, atmospheric and reactor neutrinos in the framework 
of three-neutrino mixing; in Sec. Ill we describe the data samples, the computation of 
the theoretical observables, and the statistical analysis applied in our analysis for solar 
(Sec. Ill A), atmospheric (Sec. IIIB), and reactor lllC data. Section IV is devoted to 
our results for the three-neutrino oscillation fits to solar neutrino data. Correspondingly 
in Sec. V we describe our results for atmospheric neutrino fits by themselves and also in 
combination with the reactor data. The results for the full combined five-parameter analysis 
are described in Sec. VI. Finally in Sec. Vll we summarize the work and present our 
conclusions. 

II. THREE NEUTRINO OSCILLATION PROBABILITIES 

In this section we review the theoretical calculation of the conversion probabilities for 
solar, atmospheric and reactor neutrinos in the framework of three-neutrino mixing in or- 
der to set our notation and to clarify the approximations used in the evaluation of such 
probabilities. 

In general, the determination of the oscillation probabilities both for solar neutrinos 
and atmospheric neutrinos require the solution of the Schroedinger evolution equation of 
the neutrino system in the Sun- and/or the Earth-matter background. For a three-flavour 
scenario, this equation reads 

i^ = Hz7, H = R-HjJ-Rf + V, (3) 

where R is the orthogonal matrix connecting the flavour basis and the mass basis in vacuum 
and which can be parametrized as in Eq. (1). On the other hand Hq and V are given as 

H^ = ^diag {-Aral, , 0, Ami,) , (4) 

V = diag (±V2GFNe, 0, o) , (5) 

where z7 = (z/g, z/^, z/^), Cij = cos% and Sij = sin%. The angles % can be taken without 
any loss of generality to lie in the flrst quadrant 9ij G [0,7r/2]. We have denoted by Hq 
the vacuum hamiltonian, while V describes charged-current forward interactions in matter. 
In Eq. (5), the sign + (— ) refers to neutrinos (antineutrinos), Gp is the Fermi coupling 
constant and Ne is electron number density in the Sun or the Earth. 



In writing Eq. (1) we have set all three CP violating phases to zero. Although this is, in 
general, an approximation, it holds exact for the scheme we are adopting in the description 
of solar and atmospheric data 

Am2i < Am32 ~ Amg^ , (6) 

because in this case, as we describe below, no simultaneous effect of the two mass differ- 
ences is observable in any //-appearance transition. This is the case, for instance, for the 
hierarchical scheme, mi < 1712 <^ 1713. Notice also that for transitions in vacuum the results 
obtained apply also to the inverted hierarchical case rrii > 1712 ^ ttl^. In the presence of 
matter effects the hierarchical and inverted hierarchical cases are no longer equivalent, al- 
though as discussed in Ref. [20] the difference is hardly recognizable in the current solar and 
atmospheric neutrino phenomenology. 

A. Solar Neutrinos 

For solar neutrinos we can write the survival amplitude for z/g neutrinos of energy E at 
a detector in the Earth as 

A,, = Y: A!^ Al exp[-zm,2(L - t)/2E] . (7) 

Here Afj is the amplitude of the z/g — > Ui transition {ui is the i-mass eigenstate) from the 
production point to the Sun surface, Af^ is the amplitude of the transition Ui -^ z/g from the 
Earth surface to the detector, and the propagation in vacuum from the Sun to the surface 
of the Earth is described by the exponential. L is the distance between the center of the 
Sun and the surface of the Earth, and r is the distance between the neutrino production 
point and the surface of the Sun. Using the mass hierarchy in Eq. (6) (which also implies 
that for the evolution both in the Sun and in the Earth Am32 ^ 2^/2GpNeE,ySm'^{29i3) 
the three-flavour evolution equations decouple into an effective two-flavour problem for the 
intermediate basis [31,32] 

z/g' = cos 6'i2 z/i + sin 9i2 1^2 , t'^' = - sin 612 z/i + cos 612 1^2 , (8) 

with the substitution of A^g by the effective density 

Ne' ^ Ne cos" ^13 , (9) 



while for the evolution of the third state z/,-/ = u^ there are no matter effects. Thus, the 
survival amplitude can be simplified to 

2 

Aee = cos^ ^13 ^ Af,,Ag, exp[-im^,{L - r)/2E] + sin^ ^13 exp[-iml{L - r)/2E] . (10) 

The survival probability after averaging out the interference terms due to the higher mass 
difference Am'^2 — '^"^ii is given by: 

P,3: = cos^^i3Pj:' + sm^^i3, (11) 

where we label Pjg, the corresponding two-flavour survival probability in the (Am2i,^i2) 
parameter space but with the modified matter density in Eq. (9). The expression for this 
effective two-flavour survival probability can be expressed as: 



PX' = PlPle' + P2P2e' + 2^ P^P^P^,, P^e' COS ^ • (12) 

Here Pi = |Af/jP, while Pjg' = l^fe'P and unitarity implies that P1 + P2 = 1 and Pie' + P2e' = 
1. The phase ^ is given by 

i = 1^ + ^, (13) 

where 5 contains the phases due to propagation in the Sun and in the Earth and can be safely 
neglected. In the evaluation of both Pi and P2e' the effect of coherent forward interaction 
with the Sun and Earth matter is taken into account with the effective density in Eq. (9). 
From Eq. (12) one can recover more familiar expressions for Pjg/: 

(1) For Am2i/-E~5 x 10^^^ eV, the matter effect suppresses flavour transitions both in the 
Sun and the Earth. Consequently, the probabilities Pi and P2e' are simply the projections 
of the z/e' state onto the mass eigenstates: Pi = cos^ ^12, P2e' = sin^ ^12- In this case we are 
left with the standard vacuum oscillation formula: 

p2--vac _ ^ _ gij,2(2^^2) sin2(Am2, (L - r)/AE) (14) 

which describes the oscillations on the way from the surface of the Sun to the surface of the 
Earth. 

(2) For Am\^/ E ^ 10^^^ eV, the last term in Eq. (12) vanishes and we recover the 
incoherent MSW survival probability. For Am\^/ E ~ 10^^^ — lO^"*^^ eV^, this term is zero 



because v^' adiabatically converts to V2 and Pi = 0. For Arri^i/ E^IO'^"^ eV^, both Pi and 
P2 are nonzero and the term vanishes due to averaging of cos.^. 

(3) In the intermediate range, 5 x 10^^'^ ^ Arri^i/ E ^ 10"^"' eV, adiabaticity is vio- 
lated and the cos.^ coherent term should be taken into account. The result is similar to 
vacuum oscillations but with small matter corrections. We define this case as quasi-vacuum 
oscillations [19,26,33-35,37,36]. 

In order to compute the survival probability for solar neutrinos, valid for any value of 
the neutrino mass and mixing, the full expression (11) has to be used. The results presented 
in the following sections have been obtained using the general expression for the survival 
probability in Eqs. (ll)and (12) with Pi and P2e' found by numerically solving the evolution 
equation in the Sun and the Earth matter. For Pi we use the electron number density of 
BP2000 model [38]. For P2e' we integrate numerically the evolution equation in the Earth 
matter using the Earth density profile given in the Preliminary Reference Earth Model 
(PREM) [39]. 

B. Atmospheric and Reactor Neutrinos 

For the atmospheric neutrino analysis it is a good approximation to take two of the 
neutrinos as approximately degenerate, given the hierarchy in the splittings Am^^^ and 
AttIq which is indicated by the solutions to the solar and atmospheric neutrino anomalies. 
In the Am2i -^ approximation one can rotate away the corresponding angle 612, leading 
to the following expression for the leptonic mixing matrix in vacuum [40] 



R 



/ Ci3 Si3 ^ 

-S23S13 C23 S23C13 
\ -S13C23 -S23 C23C13 J 



(15) 



^0 = ^^di^s{0,0,Aml,). (16) 

As a result the 3-neutrino propagation of atmospheric neutrinos can be well described by 
only three oscillation parameters: Aml2, ^23 and ^13. 

For ^13=0, atmospheric neutrinos involve only u^ -^ Uj. conversions, and in this case there 
are no matter effects, so that the solution of Eq. (3) is straightforward and the conversion 
probability takes the well-known vacuum form 



P,, = l-sin^(2Msin^(^^), (17) 

where L is the path-length travelled by neutrinos of energy E^. 

On the other hand, in the general case of three-neutrino scenario with ^13 7^ the 
presence of the matter potentials requires a numerical solution of the evolution equations 
in order to obtain the oscillation probabilities for atmospheric neutrinos Pap., which are 
different for neutrinos and anti-neutrinos because of the reversal of sign in Eq. (5). In our 
calculations, we will use for the matter density profile of the Earth the approximate analytic 
parametrization given in Ref . [41] of the PREM model of the Earth [39] . 

As for the CHOOZ reactor data, we need to evaluate the survival probability for v^. of 
average energy E ~ few MeV at a distance of L ~ 1 Km. For this value of energy and 
distance one can compute the survival probability neglecting Earth matter effects. In this 
case the survival probability takes the analytical form: 

pCHOOZ ^ 1 _ ^^^4 ^^^ ^i^2(2^^^) ^.^2 (^^^^ (18) 

-sin^(2^,3) (cos^^,,sin^ (^^) + ^^"' ^- ^^^^ [^^)) 
^l-sin^(2^,3)sin^(^^), 

where the second equality holds under the approximations in Eqs. (15) and (16) and is fully 

72 < 

'21 ~ 



valid for Am^, < 3 x lO^^eV^. 



III. DATA AND STATISTICAL ANALYSIS 

A. Solar Neutrinos 

In order to determine the values of neutrino masses and mixing for the oscillation solution 
of the solar neutrino problem, we have used data on the total event rates measured in the 
Chlorine experiment at Homestake [6], in the two Gallium experiments GALLEX+GNO and 
SAGE [7-9] and in the water Cerenkov detectors Kamiokande [42] and Super-Kamiokande [5] 
shown in Table I. Apart from the total event rates, Super-Kamiokande has also measured 
the dependence of the event rates during the day and during the night and the electron 
recoil energy spectrum, all measured with their recent 1117-day data sample [5]. Although, 
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as discuss in Ref. [18,19,26,43] the inclusion of Kamiokande results does not affect the shape 
of the regions, because of the much larger precision of the Super-Kamiokande data, it is 
convenient to introduce it as in this way the number of degrees of freedom (d.o.f.) for 
the rates-only-fit is 4 — 3 = 1 (instead of zero), thus allowing for the determination of a 
well-defined Xmin confidence level (CL). 

For the calculation of the theoretical expectations we use the BP98 standard solar model 
of Ref. [44]. The general expression of the expected event rate in the presence of oscillations 
in experiment i in the three-neutrino framework is given by Rj^ : 

RT = H^k fdE,Xk{E,) X [ae,.(E,)(P,,^,J (19) 

A:=l,8 



where E^, is the neutrino energy, 0^ is the total neutrino fiux and A^ is the neutrino energy 
spectrum (normalized to 1) from the solar nuclear reaction k with the normalization given 
in Ref. [44]. Here o"e,j {cx,i) is the z/g {ux, x = fi, t) interaction cross section in the SM with 
the target corresponding to experiment i. For the Chlorine and Gallium experiments we use 
improved cross sections a^^iiE) from Ref. [45]. For the Kamiokande and Super-Kamiokande 
experiment we calculate the expected signal with the corrected cross section as explained 
below. Finally {Pve^va) is the time-averaged Ve survival probability in Eq. (11). 

We have also included in the fit the experimental results from the Super-Kamiokande 
Collaboration on the day-night variation of the event rates and the recoil electron energy 
spectrum. In previous works [24,25] the data on the zenith angular dependence taken on 
5 night periods and the day averaged value, and the daily average recoil energy spectrum 
were included in order to statistically combine the information on the day-night variation 
and the energy dependence. In principle, such analysis should be taken with a grain of 
salt as these pieces of information are not fully independent; in fact, they are just different 
projections of the double differential spectrum of events as a function of time and energy 
and can be subject to possible correlations between the uncertainties in the energy and time 
dependence of the event rates, which are neglected. Here, instead, we follow the analysis of 
Ref. [26] and, in order to combine both the day-night information and the spectral data we 
use the separately measured recoil electron energy spectrum during the day and during the 
night which is free of the unknown correlated uncertainties as they correspond to different 
data samples. This will be referred in the following as the day-night spectra data which 



contains 2 x 18 data bins, including the results from the LE analysis for the 16 bins above 
6.5 MeV and the results from the SLE analysis for the two low energy bins below 6.5 MeV. 
The general expression of the expected rate in the bin j during the Day (Night) in the 
presence of oscillations, -R^^' • is similar to that in Eq. (19), with the substitution of the 
cross sections with the corresponding differential cross sections folded with the finite energy 
resolution function of the detector and integrated over the electron recoil energy interval of 
the bin, Tmin < T <T^^ 



max • 



a^ME.) = r^dT I --"^^^ dT' Res{T, T') '""'"j^;' ^'^ . (20) 



The resolution function Res{T, T') is of the form [4,46]: 

(T - T'f 



Res(T, T') = — exp 

v^(0.477T'(MeV)) 



0.44T'(MeV) 



(21) 



and we take the differential cross section daa{Ey, T')/dT' from [47]. When computing the 
spectrum during the day no Earth regeneration effect is included in the computation of 
Pue^va while during the night such effect is included as described in Sec. II. 

In the statistical treatment of all these data we perform a x^ analysis for the different 
sets of data and we define a x^ function for the set of observables X0,ratcs ^'^'^ X0,spec dn- F^^ 
the rates we follow closely the analysis of Ref. [48] with the updated uncertainties given in 
Refs. [44,45], as discussed in Ref. [24,26]. For the day-night spectra we adopt a definition 
following Ref. [24]: 

T3th nth 

X0,spec DN = Z^ Z^ [Oisp,dn pBP98 ~ « ^'^ij v'^^P"pBP98 ~ J ) \ ) 

D,Ni,j=l,18 ^i ^j 

where 

<^ij = ^ij\(^i^stat "T ^i,uncorr) "•" '^i,exp<^j,exp + (^i,cal(^j,cal (2dj 

describes the correlated and uncorrelated errors included in the Super-Kamiokande spectra 
described in Refs. [24] . Notice that in Xq ^ j-,p^ we allow for a free normalization in order 
to avoid double-counting with the data on the total event rate which is already included 
in xl j,^^gg. In the combinations of observables we define the x^ of the combination as the 
sum of the two x^'s- As discussed in Sec. II for the analysis of solar neutrino data the 
oscillation probabilities depend only on three parameters: 6'i2, 6*13 and /S.m\i. Minimizing 
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Xq obs f*^^ ^ given combination, OBS, of solar neutrino experiments as a function of the three 
neutrino oscillation parameters we determine their best fit value as well as the corresponding 
90 (95) [99] % CL allowed regions for three degrees of freedom by the condition 

x|,oBs(Am2,, ^12, ^13) - xLn,0,OBS < Ax'(CL, 3 d.o.f.) (24) 

where, for instance, Ax^(CL, 3 d.o.f.) = 6.25, 7.81, and 11.34 for CL = 90, 95, and 99% 
respectively. 

B. Atmospheric Neutrinos 

Underground experiments can record atmospheric neutrinos by direct observation of 
their charged current interaction inside the detector. These so-called contained events, can 
be further classified into fully contained events, when the charged lepton (either electron 
or muon) produced by the neutrino interaction does not escape the detector, and partially 
contained muons when the latter, produced inside, leaves the detector. For Kamiokande and 
Super-Kamiokande, the contained data sample is further divided into sub-GeV events with 
visible energy below 1.2 GeV, and multi-GeV events, with lepton energy above this cutoff. 
Sub-GeV events arise from neutrinos of several hundreds of MeV, while multi-GeV events 
are originated by neutrinos with energies of several GeV. 

Contained events have been recorded at six underground experiments, using water- 
Cerenkov detectors: Kamiokande [11], 1MB [10] and Super-Kamiokande [1,3], as well as 
iron calorimeters: Frejus [28], NUSEX [29] and Soudan2 [12]. The expected number of 
e-like and /x-like contained events, N/3 {j3 = e,fi), is given as 

Np = ntT Y. ,^ ? >^a{K c., E,)P^^-^ep{Ei) dE, dEi dc, dh (25) 

where P^I^ (-Pq^) is the v^ —^ ^p {i^a -^ ^13) conversion probability for given values of the 
neutrino energy E^, the cosine Cv of the angle between the incoming neutrino and the 
vertical direction, and the slant distance h from the production point to the sea level. In 
the SM one has P^ = 8^13 for all a, 13. In Eq. (25) rit is the number of targets, T is the 
experiment running time and <I>^ (*^a) is the fiux of atmospheric neutrinos (antineutrinos) 
of type a = e,/i, for which we will adopt the Bartol fiux [49]; Ei is the energy of the final 
charged lepton of type P = e,fi, e^{Ei) is the detection efficiency for such lepton and a^ 
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(o"^) is the neutrino- (antineutrino-) nucleon interaction cross section. Finally, Ka is the 
slant distance distribution, normalized to one [50]. For the angular distribution of events 
we integrate in the corresponding bins in q = cos 9i where 9i is the angle of the detected 
lepton, taking into account the opening angle between the neutrino and the charged lepton 
directions as determined by the kinematics of the neutrino interaction. On average the 
angle between the final-state lepton and the incoming neutrino directions ranges from 70° 
at 200 MeV to 20° at 1.5 GeV. One must also take into consideration that the neutrino 
fluxes, especially in the sub-GeV range, depend on the solar activity. In order to take this 
fact into account, we use in Eq. (25) a linear combination of atmospheric neutrino fluxes 
(^±,max g^j^j (^±,min ^^\^[Q\^ correspoud to the most active Sun (solar maximum) and quiet Sun 
(solar minimum), respectively, with different weights depending on the running period of 
each experiment [22] . The agreement of our predictions with the experimental Monte Carlo 
predictions was explicitly verified in Ref. [22]. This renders confidence in the reliability of 
our results for contained events. 

Higher energy muon neutrinos and anti-neutrinos are detected indirectly by observing 
the muons produced by charged current interactions in the vicinity of the detector: the so 
called upgoing muons. If the muon stops inside the detector, it will be called a "stopping" 
muon, while if the muon track crosses the full detector the event is classified as a "through- 
going" muon. On average stopping muons arise from neutrinos with energies around ten 
GeV, while through-going muons are originated by neutrinos with energies around hundred 
GeV. In our analysis we will consider the latest results from Super-Kamiokande [3] and 
from the MACRO [13] experiment* on upgoing muons which are presented in the form of 
measured muon fluxes. We obtain the effective muon fluxes for both stopping and through- 
going muons by convoluting the z/q, — ^ v^ transition probability (calculated as in Sec. II) 
with the corresponding muon fluxes produced by the neutrino interactions with the Earth. 
We include the muon energy loss during propagation both in the rock and in the detector 
according to Refs. [52,53] taking into account also the effective detector area for stopping 
and through-going events. Schematically 

$.(c.).,. = / r '^'ifZ'^^ ^sAE,, cME, , (26) 



*We decided not to include Baksan [51] data because they appear inconclusive. 
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where 



(P^^ i"^ 



roo fOO roo roo 

Na / dE^Q / dE^ / dX dh 
JEu Jem Jo Jo 



2^ 'jE~dir ^^ "^ "^^^ dE ^rock[Ef,o, Ef„X) [27) 

where A^^ is the Avogadro number, Efj,o is the energy of the muon produced in the neutrino 
interaction and E^ is the muon energy when entering the detector after travelhng a distance 
X in the rock. At the relevant energies the opening angle between incident neutrino and 
outgoing muon can be neglected to a very good approximation, thus we use a common label 
Ci, to characterize both directions. Here Frock{Efj,o, E^^ X) is the function which characterizes 
the energy spectrum of the muons arriving the detector. 

In Eq. (26) A^LmimCu) = As{Efj_,Cu) + At{E^,Cu) is the projected detector area for 
internal path-lengths longer than a certain Lmin- Here As and At are the corresponding 
effective areas for stopping and through-going muon trajectories. For Super-Kamiokande 
Lmin = 7 m and we compute these effective areas using the simple geometrical picture given 
in Ref. [54]. 

In contrast with Super-Kamiokande, MACRO presents its results as muon fluxes for 
E^ > 1 GeV, after correcting for detector acceptances. Therefore in this case we compute 
the expected fluxes as in Eqs. (26) and (27) but without the inclusion of the effective areas. 
In Ref. [22] we have explicitly verified that our predictions for upgoing muons agree with 
the experimental Monte Carlo predictions from Super-Kamiokande and MACRO to the 5% 
and 1% level, respectively. 

For the statistical treatment of all these data we perform a x^ analysis for different sets of 
data by computing the Xatm obs for a given combination, OBS, of experiments as a function 
of the neutrino oscillation parameters. Following closely the analysis of Refs. [22,20] we 
use the previously described contained and upgoing event numbers (instead of their ratios), 
paying attention to the correlations between the sources of errors in the muon and electron 
predictions, as well as the correlations amongst the errors of different energy data samples. 
Thus we define Xatm,OBS as 

X^..,OBS- E iNp^-Nr)-i^lA + 4H)j}-iNf^-Nr), (28) 

/.JeOBS 

with 
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^DAiTH) jj = (^DA(TH),aiA) PDA(TH),al3i^, B) a d A(T H) ,I3{.B) . (29) 

/ and J stand for any combination of experimental data sets and event-types considered, 
z.e, I = {A, a) and J = {B, 13). The latin indices A, B stand for the different experiments or 
different data samples in a given experiment. The greek indexes denote e- or /i-type events, 
z.e, a, /9 = e, /x. Here Nj^ denotes the predicted number of events (or the predicted value of 
the flux, in the case of upgoing muons) calculated as discussed above, whereas Nf"^ is the 
corresponding experimental measurement. We denote by crjjj^txH) ^^^ error matrices con- 
taining the experimental (theoretical) errors and by PDA{TH)ai3{A.,B) the matrix containing 
all the correlations between the experimental (theoretical) errors of a-like events in the A 
experiment and /3-like events in B experiment, whereas (JDA{TH),a{^) is the experimental 
(theoretical) error for the number of events in the A experiment. The dimensionality of the 
error matrix depends on the combination of experiments, OBS, included in the analysis. A 
detailed discussion of the errors and correlations used in our analysis can be found in the 
appendixes of Refs. [22] both for contained and for the upgoing muon data analysis. 

As discussed in Sec. II for the analysis of any set of atmospheric neutrino data the 
oscillation probabilities depend only on three parameters: ^13, 623 and Am'^2- Minimizing 
Xatm OBS with respcct to these three parameters we determine their best fit value as well as 
the corresponding 90 (95) [99] % confidence levels (CL) allowed regions for three degrees of 
freedom by the conditions in Eq. (24) where now Xatm obs (^"^32' ^23) ^12)- 

C. Reactor Neutrinos: CHOOZ 

The CHOOZ experiment [30] searches for disappearance of z/g produced in a power station 
with two pressurized-water nuclear reactors with a total thermal power of 8.5 GW (thermal). 
At the detector, located at L ~ 1 Km from the reactors, the z/g reaction signature is the 
delayed coincidence between the prompt e"*" signal and the signal due to the neutron capture 
in the Gd-loaded scintillator. Their measured vs. expected ratio, averaged over the neutrino 
energy spectrum is 

R = 1.01 ± 2.8 %(stat) ± 2.7 %(syst) (30) 

Thus no evidence was found for a deficit of measured vs. expected neutrino interactions 
and they derive from the data exclusion plots in the plane of the oscillation parameters 
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(Am^, sin^ 2^), in the simple two-neutrino oscillation scheme. At 90% CL they exclude the 
region given by approximately Am^ > 7-10^^ eV^ for maximum mixing, and sin^(2^) = 0.10 
for large Am^. 

In order to combine the CHOOZ results with the results from our analysis of solar and 
atmospheric neutrino data in the framework of three-neutrino mixing we have first performed 
our own analysis of the CHOOZ data. Using as experimental input their measured ratio 
(30) [30] and comparing it with the theoretical expectations we define the Xchooz function. 
As discussed in Sec. II for the analysis of the reactor data the relevant oscillation probability 
depends in general on four parameters 612, Arrigi, ^13, and Am32, but for Am2i ^ 3 x 10~^ 
eV^, even in the three-neutrino mixing scenario, it only depends on the last two. We verified 
that in this case with our Xchooz function and using the statistical criteria for two degrees 
of freedom we reproduce the excluded regions given in Ref. [30] for two-neutrino oscillations 
as can be seen in Fig. 1 where we show the excluded regions at 90, 95 and 99% CL in 
the f Am32, sin^(2^i3) j plane from our analysis of the CHOOZ data defined with 2 d.o.f. 
(Ax^jjooz ~ 4.61, 6.0, 9.21 respectively). Comparing our results with Ref. [30] we find a 
very good agreement. 

IV. THREE NEUTRINO OSCILLATION ANALYSIS OF SOLAR DATA 

As explained in Sec. II, for the mass scales involved in the explanation of the solar 
and atmospheric data the relevant parameter space for solar neutrino oscillations in the 
framework of three-neutrino mixing is a three dimensional space in the variables Am2i, ^12 
and ^13. In our choice of ordering for the neutrino masses the mass-squared difference Am'^i 
is positive and the mixing angles 612 and ^13 can vary in the interval < 612 < ^ and 
< ^13 < f . In our analysis we choose to parametrize the 612 and ^13 dependence in terms 
of the tan^ 612 and tan^ ^13 variables which span the full parameter space. 

We first present the results of the allowed regions in the three-parameter space for the 
different combination of observables. In building these regions, for a given set of observables, 
we compute for any point in the parameter space of three-neutrino oscillations the expected 
values of the observables and with those and the corresponding uncertainties we construct 
the function Xq obs(^^i' ^12 1 ^13)- We find its minimum in the full three-dimensional space 
considering both MSW and vacuum oscillations as well as the transition regime of quasi- 
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vacuum oscillations on the same footing. The allowed regions for a given CL are then 
defined as the set of points satisfying the conditions given in Eq. (24). In Figs. 2-4 we plot 
the sections of such volume in the plane (Am2i, tan^(^i2)) for different values of tan^ ^13. 

Figure 2 shows the results of the fit to the observed total rates only. For the sake of 
clarity we show the regions only at 90 and 99% CL. We find that for small tan^ ^13 ^ 0.3 
both at 90 and 99% CL, the three-dimensional allowed volume is composed of four separated 
three-dimensional regions in the MSW sector of the parameter space which we denote as 
SMA, LMA and LOW solutions analogous to the usual two-neutrino oscillation picture 
and a "tower" of regions in the vacuum oscillations (VO) sector. The global minimum 
XminQ rates = 0-62 (see Table II) used in the construction of the volumes lies in the SMA 
region and for a non-vanishing value of tan^ ^13 ~ 0.07. However, as can be seen in Fig. 5, 
this has hardly any statistical significance, as Ax^ is very mildly dependent on 613 for these 
small values of tan^ ^13. 

As seen in Fig. 2 as tan^ ^13 grows we find the following behaviours of the allowed regions 

(i) For small values of ^13 the SMA region migrates towards larger Am\^ and the LMA 
region migrates towards lower mixing angles and larger Am\^. The increase of ^13 
produces that the SMA and LMA regions merge into a unique allowed intermediate 
(INT) region which disappears at larger values of tan^ ^13 

{ii) The LOW region migrates towards the second octant and lower Am^^ and then dis- 
appears 

[ill) The 99% CL region for vacuum solution first grow for small values of ^13 and then 
become smaller as tan^ ^13 increases until it finally disappears. 

{i) and {ii) are in agreement with the results of Ref. [18] while {iii) agrees with the results 
of the second reference in [19]. 

Thus from Fig. 2 we find that as tan^ ^13 increases all the allowed regions from the fit 
to the total event rates disappear, leading to an upper bound on tan^ ^13 for any value of 
Am|;^, independently of the values taken by the other parameters in the three-neutrino 
mixing matrix. In Fig. 5.a we plot Ax|,rates(tan2 ^13) = x|,ratcs(tan^ ^13) - xLn,o,ratcs where 
X0,ratcs(^an^ 6*13) is obtained by minimizing x|,rates(^"^2i5 ^isj ^13) with respect to Am\^ and 
9i2 for fixed values of ^13 and where X^in © rates — 0-62 is the global minimum in the full 



15 



three parameter space. From the figure we can extract the upper hmit on tan^ ^13 from the 
analysis of the total event rates. The corresponding 90 and 99% CL bounds are tabulated 
in Table IV. 

Figure 3 shows the region excluded at 99% CL by the Super-Kamiokande day-night spec- 
tra data for the same tan^ ^13 values as in Fig. 2. The global minimum Xmin,&,spec dn = 28.0 
(see Table II) used in the construction of the volumes lies in the VO region and corresponds 
to tan^ 013 = 0. As seen in the figure for tan^ 613 = the excluded region overlaps with the 
lower part of the LMA allowed region (where a larger Earth regeneration effect is expected), 
covers a big fraction of the SMA region and practically the full VO region (where larger 
distortion of the energy spectrum is expected). As tan^^i3 increases the excluded region 
becomes smaller. This arises from the fact that in the survival probability (see Eq. (11)) 
the energy independent term sin^ ^13 increases and the fiat recoil electron energy spectrum 
can be more easily accommodated. For intermediate tan^ ^13 values (see third and fourth 
panels in Figs. 2 and 3 ) there is a large overlap of the excluded region with the INT region 
where the LMA and SMA merge. 

In Fig. 4 we show the results from the global fit of the full solar data-set including 
the total observed rates and the Super-Kamiokande day-night spectra data. The global 
minimum Xmm,0,giobai = 35.2 (see Table II) used in the construction of the volumes lies in 
the LMA region and corresponds to tan^ ^13 = 0. The behaviour of the regions illustrate 
the "tension" between the data on the total event rates which favour smaller ^13 values and 
the day-night spectra which allow larger values. It can also be understood as the "tension" 
between the energy dependent and constant pieces of the electron survival probability in 
Eq. (11). As a consequence of this "tension" between the two behaviours the bound on ^13 
from the global analysis, which we list in Table IV, is weaker than the one derived from 
the analysis of the event rates only. One may wonder about the meaning of the "allowed 
regions" for such large values of ^13. To clarify this point we have defined the following 
"sectors" in the (Am^i, tan^ ^12) plane 



LMA: 



1 X 10"^ < AmlJeV^ < 1 x 10"^ (31) 

1 X 10^^ < tan2^i2 < 10 
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SMA: 



1 X 10"^ < Am2^/eV2 < 1 x 10"^ (32) 

1 X 10"^ < tan^^is < 1 X 10"^ 



LOW: 



1 X 10"^ < Am2^/eV2 < 1 x 10^^ (33) 

1 X 10"^ < tan2^i2 

and we have studied the behaviour of the x© global ^^ each of these sectors. We find that in 
each of these sectors there is a local minimum around which there exists an allowed region 
if ^13 is below certain limit. As a self-consistency check we have verified that those local 
minima are well defined for any value of ^13 below the limit arising from the analysis of 
atmospheric and/or reactor data (see Sec. V). 

In Fig. 5.b we plot the values of the X0,giobai ^^ the local minima as a function of tan^ 613. 
The different curves correspond to the functions 

^X0,global(t^'^ c'i3)|scctor = Xo.globall^^'^ ^isjlscctor — Xmm,0,global 

where X© global ~ 35.2 is the global minimum in the full three-neutrino parameter space 
(see Table II) and Xo.giobailtan^ ^i3)|scctor is obtained my minimizing x|,giobai(^"^2i) ^12, ^13) 
with respect to Arrigi and 9^ in each of the above defined sectors for fixed values of ^13. In 



what follows we label as "constrained" the results of analyses when the parameters Am^ 
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and 612 are varied in a given sector while "unconstrained" refers to the case where we allow 
the variation of Arrigi and 612 in the full plane. In Fig. 5.b the curves for the constrained 
analyses are displayed up to the value of ^13 which allwos the existence of the minimum in 
the defined sectors. 

In Fig. 5.b, we also plot the function Axl^giobail^an^ ^13) = x|,giobai(tan2 ^i3)-xLn,0,giobai 
with X© global (tan^ ^13) obtained by minimizing in the full Am2i and 612 plane, i.e., for the 
unconstrained fit. This curve is simply the lower envolvent of the curves from the constrained 
analysis. From the figure we see that, unlike for the analysis of the rates only, the function 
^X© global (tan^ ^13) when the minimization is unconstrained, is not a monotonically growing 
function but presents some local maxima and minima. This behaviour is due to the fact 
that the unconstrained minimum in the (Am2i, ^12) plane moves from one sector to another: 
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• For tan^ ^13 < 0.1 the minimum lies in the LMA region. The best global fit point 
for LMA corresponds to the case tan^ ^13 = where solar and atmospheric analyses 
decouple. 

• For 0.1 < tan^ ^13 < 0.75 the minimum lies in the SMA region for which the local best 
fit point happens at tan^^i3 = 0.16. 

• At tan^ ^13 > 0.75 the minimum moves from the SMA to the INT region (which 
contains the parameter space between the three regions defined above). The preferred 
tan^ ^13 for this INT region is tan^ ^13 = 1.3. 

For larger values of tan^ ^13, a quick worsening of the x^ is produced due to the increase of 
the constant term in the survival probability. 

We can now describe more precisely the behaviour of the allowed regions shown in Fig. 4: 
for small tan^ ^13 the global fit excludes all the regions of the oscillation (Am2i, tan^ ^12) 
parameter plane where the energy dependent piece of the survival probability in Eq. (11), 
cos^ OisP^,"^,, is either too little to account for the observed event rates or too large to account 
for the flat spectrum. As ^13 increases the constant sin^ ^13 piece in the survival probability 
increases and as a result the fit to the fiat spectrum becomes good in the full (Am25^, tan^ 612) 
plane while the fit to the event rates becomes better in those regions where the energy 
dependence of the Pjg/ piece is stronger. In this way for tan^ ^13 ^ 0.4 the local minimum 
in the (Am2n tan^ 612) plane for fixed tan^ ^13 moves from the LMA region first to the SMA 
and finally to the INT region where for low energy neutrinos the energy dependence of Pjg/ 
is stronger. As tan^ ^13 increases to much larger values the survival probability becomes 
basically energy independent and the fit to the event rates becomes too bad in the full 
parameter space. 

Let us finally comment on the statistical meaning of the allowed regions. Notice that 
following the standard procedure the allowed regions shown in Figs. 2-4 have been defined in 
terms of shifts of the x^ function for those observables with respect to the global mininium. 
Defined this way, the size of a region depends on the relative quality of its local minimum 
with respect to the global minimum but from the size of the region we cannot infer the 
actual absolute quality of the description in each region. That is given by the value of the 
X© global function at the local minimum (which for this case we show in Fig. 5.b). From 
this analysis we see that for small tan^ 6*13 the values of Xo.giobai ^^ ^^^ local minimum in 



the different regions are not so different. For instance, for tan^ ^13 < 0.2 we find that the 
goodness of the fit (GOF) for the different solutions is 55% for the LMA and 37% for the 
SMA and LOW solutions^. Thus our conclusion is that from the statistical point of view 
for small tan^ ^13 all solutions are acceptable since they all provide a reasonable GOF to the 
full data set. Although LMA solution seem slightly favoured over SMA and LOW solution 
these last two solutions cannot be ruled out at any reasonable CL. 

V. THREE NEUTRINO OSCILLATION ANALYSIS OF ATMOSPHERIC AND 

REACTOR DATA 

A. Atmospheric Neutrino Fit 

In our statistical analysis of the atmospheric neutrino events we use the following 
data: (i) unbinned contained e-like and /i-like event rates from Frejus [28], 1MB [10], Nu- 
sex [29], Kamiokande sub-GeV [11] and Soudan2 [12]; (ii) e-like and /x-like data samples 
of Kamiokande multi-GeV [11] and Super-Kamiokande sub- and multi-GeV [3], each given 
as a 5-bin zenith-angle distribution^; (iii) Super-Kamiokande upgoing muon data including 
the stopping (5 bins in zenith-angle) and through-going (10 angular bins) muon fluxes; (iv) 
MACRO [13] upgoing muons samples, with 10 angular bins. 

In order to study the results for the different types of atmospheric neutrino data we have 
deflned the following combinations of data sets: 

• FINKS The e-like and /x-like event rates from the flve experiments Frejus, 1MB, Nusex, 
Kamiokande sub-GeV and Soudan-2. It contains 10 data points. 

• CONT-UNBIN The rates in FINKS together with Kamiokande multi-GeV and Super- 
Kamiokande sub and multi-GeV e-like and /x-like event rates without including the 
angular information, which accounts for a total of 16 data points. 



TThe small differences in the GOF with the results in Ref. [26] are due to the effect of the additional 
^13 parameter 

■•■Note that for convenience and maximal statistical significance we prefer to bin the Super- 
Kamiokande contained event data in 5, instead of 10 bins. 
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• CONT-BIN The rates in FINKS together with Kamiokande multi-GeV and Super- 
Kamiokande sub and muhi-GeV e-hke and /i-hke event rates including the angular 
information. It contains 40 data points. 

• UP-/i Muon fluxes for stopping and through-going muons at Super-Kamiokande and 
MACRO which correspond to 25 data points. 

• SK The angular distribution of e-like and /x-event rates and upgoing muon fluxes 
measured at Super-Kamiokande. It contains 35 data points. 

• ALL-ATM The full data sample of atmospheric neutrino data which corresponds to 
65 points. 

The first result of our analysis refers to the no-oscillation hypothesis. As can be seen from 
the fifth column of Table III, the x^ values in the absence of new physics - as obtained with 
our prescriptions for different combinations of atmospheric data sets - clearly show that the 
data are totally inconsistent with the SM hypothesis. In fact, the global analysis, which refers 
to the full combination ALL-ATM, gives Xsm all-atm ~ 191.7/(65 d.o.f.) corresponding to 
a probability ^ 10^^^. This result is rather insensitive to the inclusion of the CHOOZ reactor 
data, as can be seen by comparing the values of x^ given in the last two lines of Table III. 
This indicates that the Standard Model can be safely ruled out. In contrast, the x^ for 
the global analysis decreases to xLin = 61.7/(62 d.o.f.) [Xmin = 62.5/(63 d.o.f.) including 
CHOOZ], acceptable at the 51% when oscillations are assumed. 

Table III also gives the minimum x^ values and the resulting best fit points for the various 
combinations of data sets considered. Note that for FINKS, CONT-UNBIN and UP-/i 
combinations the best-fit point is characterized by a rather large value of tan^ ^13, while all 
the other data set favour a value very close to 0. The corresponding allowed regions for the 
(tan^ 623, Am'^2) oscillation parameters at 90, 95 and 99% CL for the different combinations 
are depicted in the set of Figures from Fig. 6-Fig. 13. In all these figures the upper-left 
panel, tan^ ^13 = 0, corresponds to pure u^ -^ v^ oscillations, and one can note the exact 
symmetry of the contour regions under the transformation ^23 -^ 7r/4 — ^23- This symmetry 
follows from the fact that in the pure v^ -^ v^ channel matter effects cancel out and the 
oscillation probability depends on ^23 only through the double-valued function sin^(2^23)- 
For non-vanishing values of ^13 this symmetry breaks due to the three-neutrino mixing 
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structure even if matter effects are neglected. With our sign assignment we find that in 
most cases for non-zero values of ^13 the allowed regions become larger in the second octant 

of 6*23. 

In Figs. 6 and 7 we present the allowed regions in (tan^^23, ^^^32) fo^^ different values 
of tan^ ^13, for the FINKS and CONT-UNBIN data respectively. It is evident that, despite 
of the large statistics provided by Super-Kamiokande data, it is not possible from the 
information on the total event rates only, without including the angular dependence, to 
place any upper bound on Am32, and even the lower bound Am32 > 2 x 10~^ eV^ is rather 
weak. The CONT-UNBIN data also do not provide a relevant constraint on tan^ ^13. 

Conversely, in Fig. 8 we display the allowed regions in (tan^^23j Am32) for different 
values of sin^^i3, for the combination of CONT-BIN events, including also the information 
on the angular distributions. Note that, as expected, the upper bound on ^171^2 is now 
rather strong (better than 10"^ eV^) as a consequence of the fact that no suppression for 
downgoing u^ neutrinos is observed. This imposes a lower bound on the neutrino oscillation 
length and consequently an upper bound on the mass difference. However contained events 
alone still allow values of ATn'^2 < 10^^ eV^. Note also that the allowed region is still rather 
large for tan^ ^13 ^ 0.7 and at 99% CL it only disappears for tan^ ^13 ~ 2.4. 

In order to illustrate the main effect of adding the angle ^13 in the description of the 
angular distribution of contained atmospheric neutrino events we show in Fig. 9 the zenith- 
angle distributions for the Super-Kamiokande e-like (left panels) and /x-like (right-panels) 
contained events, both in the sub-GeV (upper panels) and multi-GeV (lower panels) energy 
range. The thick solid line is the expected distribution in the absence of oscillation (SM 
hypothesis), while the thin full line represents the prediction for the overall best-fit point 
of the full atmospheric data set (ALL-ATM) (see also Table III) which occurs at a small 
tan2^i3 = 0.026 (with Aml^ = 3.3 x 10"^ eV^, tan^ ^23 = 1-63). The dashed and dotted 
histograms correspond to the distributions with increasing value of tan^ ^13 = 0.33, 0.54 
which are the maximum allowed values at 90 and 99% CL from the analysis of all atmospheric 
data. For each such tan^ ^13 we choose Attl'^^ s-'^d tan^ ^32 so as to minimize the x^- Clearly 
the oscillation description is excellent as long as the oscillation is mainly in the z/^ -^ Vr 
channel (small ^13). This is simply understood since, from the left panels, it is clear that 
e-like events are well accounted for within the no-oscillation hypothesis. From the figure we 
see that increasing tan^ ^13 leads to an increase in all the contained event rates. This is due 
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to the fact that an increasing fraction of z/^ now oscillates as z/^ -^ Ve (also z/g's oscillate as 
^e ~^ ^11 but since the z/g fluxes are smaller this effect is relatively less important) spoiling 
the good description of the e-type data, especially for upgoing multi-GeV electron events. 
For multi-GeV events all the curves coincide with the SM one for downgoing neutrinos which 
did not have the time to oscillate. For sub-GeV this effect is lost due to the large opening 
angle between the neutrino and the detected lepton. We also see that for multi-GeV electron 
neutrinos the effect of ^13 is larger close to the vertical (cos^ = — 1) where the expected 
ratio of fluxes in the SM i?(z/^/z/e) is larger. Conversely the relative effect of ^13 for z/^ is 
larger close to the horizontal direction, cos^ = 0. 

Now we move to upgoing muon events. In Fig. 10 we show the allowed regions in 
(tan^^235 ^^^32) for different tan^ ^13 values, for the combination UP-/i which contains 
upgoing-muon events from Super-Kamiokande (through-going and stopping) and MACRO 
(through-going only). This plot is complementary to Fig. 8 (corresponding to the CONT- 
BIN combination), in the sense that the data in CONT-BIN and UP-/i combinations are 
completely disjoint. In contrast to the CONT-BIN case, now the upper bound on Am|2 is 
much weaker (^ 3 x 10~^ eV^), while the lower bound is now stronger. Again, no relevant 
bound can be put on ^13 from the analysis of upgoing events alone. 

The angular distribution for the upward-going muon fluxes for increasing values of ^13 
is presented in Fig. 11. The thick solid line is the expected distribution in the absence of 
oscillations (SM hypothesis), while the thin full line represents the prediction for the overall 
best-fit point of all atmospheric data (ALL-ATM). As in Fig. 9 the dashed and dotted 
histograms correspond to the distributions with increasing value of tan^ ^13 = 0.33, 0.54 
(maximum acceptable values at 90 and 99% CL from the analysis of ALL- ATM data). 
From the figure we see that the effect of adding a large ^13 in the expected upward muon 
fluxes is not very significant. For stopping muons the effect is larger for neutrinos arriving 
horizontally. This is due, as for the case of multi-GeV muons, to the larger i?(z/e/z/^) SM 
flux ratio in this direction which implies a larger relative contribution from z/g oscillating 
to u^. This feature is lost in the case of through-going muons because this effect is partly 
compensated by the matter effects and also by the increase of tan^ ^23- 

In order to perform a separate critical analysis of the implications of all Super- 
Kamiokande data by themselves we show in Fig. 12 the allowed regions in (tan^ ^231 ^^^32) 
for different tan^ ^13 values, for the combination of SK data (contained and upgoing). Due 
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to the large statistics provided by this experiment, Am'^2 is strongly bounded both from 
above and from below. Moreover no region of parameter space is allowed (even at 99% CL) 
for tan^ ^13 > 0.7. It is also interesting to notice that (unlike in the ALL- ATM combination 
discussed later) as tan^ ^13 increases the allowed region in the second octant of 623 becomes 
smaller and finally disappears. This behaviour is driven by the SK contained event data 
which favours the first 623 octant as can be seen in the corresponding line in Table IIL 

Finally we discuss the results from the combined global analysis of all the atmospheric 
neutrino data (ALL- ATM). The allowed range of parameters are shown in Figs. 13 and 16. 
In Figs. 13 we show the global (tan^^23, ^^^32) allowed regions, for different values of 
tan^ ^13 while in Fig. 16 we show the corresponding projection of the three-dimensional 
parameter space in the (tan^^23, tan^ ^13) plane, for different values of Am|2- Although 
Fig. 13 shows no qualitative difference with respect to the allowed regions from the analysis 
of SK data alone displayed in Fig. 12, we find that the inclusion of the other experimental 
results still results into an slightly tighter restriction on the allowed parameter space. In this 
way, for instance, we find that the allowed region (at 99% CL) disappears for tan^ ^13 ~ 0.6. 
Moreover, as mentioned before comparing Fig. 12 with 13 we see a slight trend in Fig. 12 
towards the tan^ 623 < 1 octant, while tan^ 623 > 1 is favoured in the other case. 

All these features can be more quantitatively observed in Figs. 14 and 15 where we show 
the dependence of the Axltm function on tan^ ^13 and on ^171^2, for the different combina- 
tion of atmospheric neutrino events. In these plots all the neutrino oscillation parameters 
which are not displayed have been "integrated out", or, equivalently, the displayed Ax^ is 
minimized with respect to all the non-displayed variables. From the left panels of Figs. 14 
and 15 we can read the upper bound on tan^ ^13 that can be extracted from the analysis 
of the different samples of atmospheric data alone regardless of the values of the other pa- 
rameters of the three-neutrino mixing matrix. The corresponding 90 and 99% CL bounds 
are listed in Table IV. Conversely from the right panels of Figs. 14 and 15 we extract the 
allowed value of A'm'^2 by the different combination irrespective of the of the values of the 
other parameters of the three-neutrino mixing matrix. 

In conclusion we see that the analysis of the full atmospheric neutrino data in the frame- 
work of three-neutrino oscillations clearly favours the z/^ -^ u^- oscillation hypothesis. As a 
matter of fact the best fit corresponds to a small value of ^13 = 9°. But it still allows for a 
non- negligible z/^ -^ Ve component. More quantitatively we find that the following ranges of 
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parameters are allowed at 90 [99] % CL from this analysis 

[1.25 X 10"^] 1.6 X 10"^ < Amla/eV^ < 6 x 10"^ [8 x 10"^] (34) 

[0.37] 0.43 < tan2^23 < 4.2 [6.2] 
tan2^i3 < 0.34 [0.57] 

One must take into account that these ranges are strongly correlated as illustrated in Figs. 13 
and 16. 

B. Fit to Atmospheric and CHOOZ data 

We now describe the effect of including the CHOOZ reactor data together with the 
atmospheric data samples in a combined 3-neutrino x^ analysis. The results of this analysis 
are summarized in Fig. 15, Figs. 17-19, and in Tables. IV and III. In this analysis we will 
assume that ^m^^ ^ 3 x 10^^ eV^ and work under the approximations in Eqs. (15) and (16). 

As discussed in Sec. Ill C the negative results of the CHOOZ reactor experiment 
strongly disfavour the region of parameters with Am32 ^ 10~^ eV^ and sin^(2^i3) ^ 0.10 
(0.026^ tan^ ^13^38). However for smaller values of Am32 the CHOOZ result leads to much 
weaker bounds on the ^13 mixing angle. To illustrate this point we show in Fig. 17, the 
allowed regions from the combination of the CONT-BIN events with the CHOOZ data. We 
see in Fig. 17, which should be compared with Fig. 8, that, as soon as tan^^i3 deviates 
from zero the ^m\2 > 10~^ eV^ region is ruled out. However there is still a region in the 
parameter space which survives (at 99% CL) up to tan^ ^13 ^ 0.66. Thus, even with the 
large statistics provided by the Super-Kamiokande data and including the CHOOZ result, 
it is not possible to constrain ^13 using only the data on contained events. 

The situation is changed once the upgoing muon events are included in the analysis. As 
shown in Figs. 10 and Fig. 13, the upgoing muon data disfavours the low mass Aml2 < 10^^ 
eV^. As a result the full 99% CL allowed parameter regions from the global analysis of the 
atmospheric data shown in Fig. 13 lies in the mass range where the CHOOZ experiment 
should have observed oscillations for sizeable ^13 values. This results into the shift of the 
global minimum from the combined atmospheric plus CHOOZ data to ^13 = 0. Thus adding 
the reactor data (Fig. 18) has as main effect effect the strong improvement of the tan^ ^13 
limit, as seen from the left panel in Fig. 15 and by comparing the allowed ranges in Fig. 16 
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and Fig. 19. From these figures we see that the higher Am^2 the more the CHOOZ data 
restricts ^13. In this way the allowed range of ^13 at 99% CL is reduced by a factor ~ 4, 10, 
15, 25 for Am32 = 1.5,3.0,4.5 and 6xl0~^ eV^ respectively by the inclusion of the CHOOZ 
data. In contrast, the allowed range of A?72|2 is only mildly restricted when combining the 
full atmospheric data with the reactor result as displayed in the right panel of Fig. 15. 

We get as final results from the atmospheric and reactor neutrino data analysis the 
following allowed ranges of parameters at 90 [99] % CL 

[1.2 X 10"^] 1.6 X 10"^ < Aml^/eV^K 5.4 x 10"^ [6.6 x 10"^] (35) 

[0.36] 0.43 < tan^ 623 < 2.7 [3.3] 

tan2 0i3 < 0.043 [0.08] 

Comparing with the allowed ranges shown in Eq. (35) we see that the main effect of the 
addition of the CHOOZ results to the atmospheric neutrino data analysis is the stronger 
constraint on the allowed value of the angle ^13. This illustrates again the fact that the full 
allowed region of Arrigg from the atmospheric neutrino analysis lies in the sensitivity range 
for the CHOOZ reactor. However it only holds once upgoing niuons are included in the fit. 
Including the CHOOZ results has very little effect on the allowed range of Am|2 while it 
results into a tighter constraint on tan^ ^23 on the second octant. This last effect arises from 
the fact that in order to have tan^ 623 far into the second octant one requires large ^13 values 
which are forbidden by the CHOOZ data. 

VI. COMBINED SOLAR, ATMOSPHERIC AND REACTOR ANALYSIS 

In this section we describe the results of the combined analysis of the solar and atmo- 
spheric neutrino data by themselves and also in combination with the reactor results. In 
order to perform such an analysis we have added the x^ functions for each data set. In this 
way we define 

Xatm + solar(^"^21> ^"^32> ^12, ^23, ^13) = X0,global(^"^21, ^12, ^13) + (36) 

Xaim,ALL(^"^32) ^23, ^13) 
Xatm + solar + CHOOz(^"^21. ^"^32' ^12, ^23, ^13) = X0,global(^"^21) ^12, ^13) + 

Xatm,ALL(^"^32>^23,^13) + 

+XcHooz(^"^2i, ^12, Am^2, ^13). 
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Notice that the x^ functions defined in this way depend on 5 parameters. Therefore the 
allowed parameter space at a given CL is a 5-dimensional volume defined by the corre- 
sponding conditions on A^^ for 5 d.o.f., A^^ < 9.24 (11.07) [15.09] at 90 (95) [99] % CL. 
The results of the global combined analysis are summarized in Fig. 20 and Fig. 21 and in 
Tables V and VI. 

In Fig. 20 we plot the behaviour of the Ax^ functions defined above with respect to 
At7i|2, tan^ ^23 and tan^^i3. In the upper panels we show the behaviour when including 
in the analysis only the solar and atmospheric neutrino data while in the lower panels the 
CHOOZ reactor constraint is added. 

In constructing these Ax^ functions we have subtracted the minimum in the 5-parameter 
space and we have minimized with respect to the "solar" parameters, Arrigi and tan^ 612, as 
well as the other two undisplayed atmospheric parameters. Both the 5-parameter minimum 
and the solar minimization have been performed either in the full (Am2i, tan^ 612) parameter 
plane (which we previously defined as unconstrained) or restricting the solar parameters to 
lie in the LMA or SMA sectors (as defined in Eqs. (32) and (33), see also the discussion 
below on the solar parameter space). The corresponding values of Xmin ^^^ ^^e position 
of the minima in the 5-dimensional space for each case are given in Table V. We have 
verified that for all curves in Fig. 20 the minimization in the solar parameters always occurs 
at Am2i < 10^'^ eV^ which ensures the validity of the approximations in Eqs. (15) and (16). 

Let us first discuss the effect of the combination of solar, atmospheric and reactor data 
on the common parameter to the full analysis, ^13. In Fig. 20. a we display the dependence 
of Axit^ _|_ solar with tau^ ^13 ouce we have minimized in the other four parameter. As can be 
seen from the figure A^^^-j^j ^ g^j^j, is sensitive to the particular solution of the solar neutrino 
problem LMA, SMA or LOW. As we discussed in Sec. IV tan^ ^13 ~ 0.1 is the turning point 
where the minimum for the solar neutrino analysis switches from LMA to SMA. This change 
produces the features of the unconstrained curve. Notice that for tan^ ^13 ^0.1 the uncon- 
strained curve does not coincide with the SMA one because the 5-dimensional minimum 
subtracted is different. But, as expected, they are roughly paralell. The corresponding 90 
and 99% limits can be read in Table VI. Correspondingly Fig. 20. d shows the results once 
the CHOOZ data is added. We see that the inclusion of CHOOZ produces a tighter restric- 
tion in tan^^i3, now bounded to lie below 0.1 at 99% CL. As a result the unconstrained 
minimum is always in the LMA region. On the other hand SMA prefers finite ^13 leading 
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to a less restrictive limit on ^13 for the global fit. The allowed 90 and 99% CL ranges are 
shown in Table VI. 

Adding the solar data to the reactor and atmospheric results also affects the allowed 
ranges on the "non-solar" parameters ^23 and Aml2- This effect is illustrated in panels 
(b), (c), (e) and (f) of Fig. 20. Figure 20. b shows the dependence of Ax^^.^^^ _^ g^i^^j, with 
tan^ 623 after minimizing in the other four parameters. This figure illustrates that the allowed 
tan^ ^23 range is sensitive to the solar solution (LMA or SMA) only far in the second octant 
(tan^ ^23 > 2). As discussed in Sec. V such region of ^23 is favoured for larger values of ^13. 
For such ^13 values Axq global depends on the solar solution. Hence the behaviour observed. 
As seen in Fig. 20. e once we include the CHOOZ results, the sensitivity on the solar solution 
disappears because now we only have very small ^13 values. Comparing Fig. 20. b with 
Fig. 20. e one sees that the inclusion of CHOOZ data results in a tighter constraint on ^23 
mainly in the second octant irrespective of the allowed region for the solar parameters. On 
the other hand comparing Fig. 20. c with 20. f one sees that the lower limit on Am^2 is rather 
insensitive to the type of solar neutrino solution while the upper limit gets reduced and also 
becomes independent on the solar solution. 

Let us now discuss how the allowed range of solar parameters, Am^i and 612, is affected 
by the inclusion of the atmospheric and reactor data in the analysis. In order to illustrate 
this point we display in Fig. 21 the allowed regions in Amgi and 612 at 90 and 99% CL (for 
5 degrees of freedom) after minimizing the function Xatm + solar +chooz with respect to the 
other three parameters. We show the figure with two different statistical criteria. 

In Fig. 2 La we show the allowed regions for the solar parameters for the unconstrained 
analysis. The regions are all defined with respect to the global 5-dimensional minimum 
(which occurs in the LMA region) given in Table V. This is the same criterion used to 
define the regions of the 3-dimensional solar analysis in Figs. 2-4. However here, instead of 
showing the section for a fixed value of ^13 this parameter has also been minimized in order 
to obtain the maximum allowed regions in the (Atti^i, tan^i2) plane. As seen from the figure 
the allowed parameter space still consists of three well defined disjoint regions. This is due 
to the smallness of ^13. The main difference with the first panel in Fig. 4 arises from the 
fact that the allowed value of Ax^ at a given CL used in the definition of the regions is now 
different due to the additional degrees of freedom. Also, we see that the LMA region is cut 
for Amli ^ 7 X 10~'^ as a consequence of the inclusion of the CHOOZ data (see Eq.(19)). 
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In Fig. 21. b each of the three regions: LMA, SMA and LOW are defined with re- 
spect to their local minimum, also given in Table V. These are the parameter ranges in 
(Am^i, tan^^i2) allowed by the combined analysis, but constraining the solar analysis to a 
given sector. Notice also that since the global minimum is in LMA, this region is the same 
in Fig. 21. a and Fig. 21. b, while the SMA and LOW regions are different. 

Comparing for example the first panel of Fig. 4 with Fig. 21. a one sees how the exact 
size of the allowed solar regions at a given CL depends on the presence of additional degrees 
of freedom. This again illustrates how, as discussed in Sec. IV, we cannot infer the actual 
quality of the description in a region from its size only. 

Summarizing, our final results from the joint solar, atmospheric, and reactor neutrino 
data analysis lead to the following allowed ranges of parameters at 90 [99] % CL 

[1.1 X 10"^] 1.4 X 10"^ < ^ml^/eV^K 6.1 x lO"''^ [7.3 x lO"''^] (37) 

[0.33] 0.39 < tan2^23 < 3.1 [3.8] 

tan2 6'i3 < 0.055 [0.085] (Unconstrained or LMA) 

tan2^i3 < 0.075 [0.135] (SMA) 

On the other hand the allowed values of ^m\^ and 9i2 can be inferred from Fig. 21. a or 
Fig. 21. b for unconstrained and constrained cases respectively. 

VII. SUMMARY AND CONCLUSIONS 

In this article we have performed a three-fiavour analysis of the full atmospheric, solar 
and reactor neutrino data. The analysis contains five parameters: two mass differences, 
Am32 and Am2i, and three mixing, ^12, ^23, and ^13. Under the assumption of mass hierarchy 
in neutrino masses the solar neutrino observables depend on three of these parameters, 
which we chose to be Am'^i, ^12 and ^13 while the atmospheric neutrino event rates depend 
on Am32, ^23 and ^13. The survival probability for reactor neutrinos depends in general 
on four parameters 612, Am2i, ^13, and Am32, but for Am2i ^ 3 x 10^^ eV^ it effectively 
depends only on ^13 and Am|2- Thus we have that in the hierarchical approximation the 
only parameter common to the three data sets is ^13. 

First we have performed independent analyses of the solar neutrino data and of the 
atmospheric (and atmospheric plus reactor) neutrino data in the respective 3-dimensional 
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parameter spaces. In our solar data analysis we have studied the dependence of the solutions 
to the solar neutrino problem on the ^13 parameter. We find that the most favourable 
scenario for solar neutrino oscillations is the simplest two-neutrino mixing case, ^13 = and 
that for large enough ^13 angles there is no allowed solution to the solar neutrino problem. 
As a result we derive an upper bound on tan^ ^13 < 2.4 at 90% CL. 

Our results for the three-dimensional analysis of the full atmospheric neutrino data in 
the framework of three-neutrino oscillations show that the most favourable scenario is the 
^'^ ~^ ^T oscillation hypothesis with the best fit corresponding to a very small value of 
^13 = 9°. However a non-neglegeable z/^ -^ Ve component is still allowed. The corresponding 
upper limit on ^13 is tan^ ^13 < 0.34 at 90% CL. From this study we have also derived the 
allowed Amgg and tan^ ^32 ranges in the framework of three neutrino mixing. 

We have also studied the effect of combining the CHOOZ and the atmospheric neutrino 
data in a common three-parameter analysis. Our results show that the main effect of 
the addition of the CHOOZ data to the atmospheric neutrino analysis is to strengthen the 
constraint on the allowed value of the angle ^13 which leads to the upper limit tan^ ^13 < 0.043 
at 90% CL. This is due to the fact that the full allowed region of Am32 from the atmospheric 
neutrino analysis lies in the sensitivity range for the CHOOZ reactor. Including the CHOOZ 
results has very little effect on the allowed range of Amig, while it results into a tighter 
constrain on tan^ ^23 on the second octant. 

Finally we have performed the combined five-dimensional analysis of the solar and atmo- 
spheric neutrino data and also in combination with the reactor results and we have derived 
the allowed range of the five parameters. We have also studied how these ranges depend on 
the particular solution region for the solar neutrino deficit. 

In conclusion we see that from our statistical analysis of the solar data it emerges that 
the status of the large mixing-type solutions has been further improved with respect to 
the previous Super-Kamiokande data sample, due mainly to the substantially fiatter recoil 
electron energy spectrum. In contrast, there has been no fundamental change, other than 
further improvement due to statistics, on the status of the atmospheric data. For the latter 
the oscillation picture clearly favours large mixing, while for the solar case the preference 
is still not overwhelming. Both solar and atmospheric data favour small values of the 
additional ^13 mixing and this behaviour is strengthened by the inclusion of the reactor 
limit. Nevertheless the prospects that both solar and atmospheric data select large lepton 
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mixing seems in puzzling contrast with the observed structure of quark mixing. 
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TABLES 



Experiment 


Rate 


Ref. 


Units 


RfP^^ 


Homestake 


2.56 ±0.23 


[6] 


SNU 


7.8 ±1.1 


GALLEX + SAGE+ GNO 


74.7 ±5.2 


[8,7] 


SNU 


130 ±7 


Kamiokande 


2.80 ± 0.38 


[42] 


10^ cm-2 s-i 


5.2 ±0.9 


Super-Kamiokande 


2.40 ± 0.08 


[5] 


10^ cm-2 s-i 


5.2 ±0.9 



TABLE L Measured rates for solar neutrinos at by Chlorine, Gallium, Kamiokande and Su- 
per-Kamiokande experiments. 



Data sets 


d.o.f. 


tan2 6'i3 tan2 6'i2 Am^i [eV^] 


v2 v2 
ASM Amin 


Rates 

SpecDAT 

Global 


4-3 
35-3 
39-3 


0.07 1.0 X 10-3 8.2 X 10-6 
0.0 0.17(5.9) 4.6 X 10-^° 
0.0 0.37 3.3 X 10-5 


62.2 0.62 
30.4 28.0 
92.6 35.2 



TABLE IL Minimum x^ values and best-fit points for various sets of solar neutrino data. 



Data sets 


d.o.f. 


tan^ 6*13 


tan^ 023 


Ami, [eV^] 


XSM 


Amin 


FINKS 


10-3 


1. 


> 100 


0.9 X 10-3 


29.9 


14.8 


CONT-UNBIN 


16-3 


0.33 


> 100 


4.2 X 10-3 


58.1 


18.4 


CONT-BIN 


40-3 


0.03 


1.6 


2.5 X 10-3 


115.1 


41.3 


SK CONT 


20-3 


0.03 


0.89 


2.3 X 10-3 


80.7 


14.5 


UP-^ 


25-3 


0.23 


3.3 


7.1 X 10-3 


50.3 


25.1 


SK 


35-3 


0.005 


1.3 


2.8 X 10-3 


131.4 


27.8 


ALL-ATM 


65-3 


0.03 


1.6 


3.3 X 10-3 


191.7 


61.7 


ALL-ATM ± CHOOZ 


66-3 


0.005 


1.4 


3.1 X 10-3 


191.8 


62.5 



TABLE III. Minimum x^ values and best-fit points for various sets of atmosperic neutrino data 
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tan^ ^13 


6*13 deg 




data set 


mill 


limit 90% (99%) 


limit 90% (99%) 


Solar 


Rates 


0.22 


0.8 (2.5) 


42° (58°) 




Global 


0.0 


2.4 (3.5) 


57° (62°) 




FINKS 


L 


16.4 (42.7) 


76° (81°) 




CONT-UNBIN 


0.33 


10.0 (13.6) 


72° (75°) 




CONT-BIN 


0.03 


0.97 (2.34) 


45° (57°) 




UP-/i 


0.23 


0.55 (7.71) 


37° (70°) 




SK 


0.005 


0.33 (0.68) 


30° (40°) 




ALL-ATM 


0.03 


0.34 (0.57) 


30° (37°) 


ALL-ATM + Chooz 


0.005 


0.043 (0.08) 


12° (16°) 



TABLE IV. Upper bounds on ^13 at 90 and 99% CL from the analysis of the different smaples 
and combinations of atmosperic neutrino data. 



data set 


Solar region 


tan^ 6*13 


tan^ 6*12 


Amii [eV^] 


tan^ 6*23 


Ami2 [eV2] 


Amin 


Solar + ALL-ATM. 


Unconst=LMA 


0.015 


0.35 


3.3 X 10-5 


1.6 


3.3 X 10-3 


97.1 




LOW 


0.025 


0.70 


9.2 X 10-^ 


1.6 


3.3 X 10-3 


100.8 




SMA 


0.03 


5.8 X 10-4 


5.6 X 10^6 


1.7 


3.3 X 10-3 


101.5 


Solar + ALL-ATM. 


Unconst=LMA 


0.005 


0.36 


3.3 X 10-5 


1.4 


3.1 X 10-3 


97.7 


+CHOOZ 


LOW 


0.005 


0.58 


9.6 X 10"^ 


1.4 


3.1 X 10-3 


101.4 




SMA 


0.005 


6.8 X 10-4 


5.1 X 10-^ 


1.4 


3.1 X 10-3 


103.1 



TABLE V. Minimum x^ values and best-fit points for various sets the combined analysis of 
atmospheric, solar and reactor neutrino data. 
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limit 90% (99%) 



data set 



Solar region 



tim^e 



13 



tan^ 6*23 



[1.35, 6.5]([1. 1,8.0]) X 10-3 
1.35, 6.5]([1.07, 8.0]) x lO^^ 
[1.35, 7.5]([1. 1,9.7]) X 10-3 



Solar + ALL-ATM. 



Unconstrained 
LMA 
SMA 



< 0.32(< 0.52) 

< 0.20(< 0.39) 

< 0.46(< 0.75) 



[0.39, 4.3] ([0.33, 6.5]) 
[0.39, 3.6]([0.33, 5.2]) 
[0.40, 6.0] ([0.34, 8.3]) 



Solar + ALL-ATM 
+ CHOOZ 



Unconst=LMA 
SMA 



< 0.055(< 0.085) 

< 0.075(< 0.135) 



[0.39, 3.0] ([0.33, 3.7]) 
[0.39, 3.1]([0.33, 3.8]) 



[1.4, 6.1]([1. 1,7.3]) X 10-3 
[1.3,6.1]([1.0,7.3]) X 10-3 



TABLE VI. Allowed ranges on tan^e^is, tan^ 6*23 and Am^a [eV^] at 90 and 99% CL from the 



combined analysis. 
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FIG. 1. Excluded region in Arrigj and sin (2^13) from the non observation of oscillations by the 
CHOOZ reactor experiment. The curves represent the 90, 95 and 99% CL excluded region defined 
with 2 d.o.f. for comparison with the CHOOZ published results. 
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FIG. 2. Allowed three-neutrino oscillation regions in Am^i and tan^ 9i2 from the measure- 
ments of the total event rates at Chlorine, Gallium, Kamiokande and Super-Kamiokande (1117-day 
data sample) experiments. The different panels represent the allowed regions at 99% (darker) 
and 90% CL (lighter) obtained as sections for fixed values of the mixing angle tan^ ^13 of the 
three-dimensional volume defined by x^ — Xmm=6. 25 (90%), 11.36 (99%). The best-fit point is 
denoted as a star. 
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FIG. 3. Excluded three-neutrino solar oscillation regions at 99% CL in Atti^i and tan^ 612 from 
the measurement of the day-night spectra data by Super-Kamiokande (1117-day data sample). 
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FIG. 4. Allowed three-neutrino solar oscillation regions in Am2i and tan^ ^12 from the global 
analysis of solar neutrino data. The best-fit point is denoted as a star. 
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FIG. 5. Ax^ as a function of tan^ ^13 from the three-neutrino analysis of the solar data. The 
left panel corresponds the analysis of total rates only and the right panel to the global analysis. 
The dotted horizontal lines correspond to the 90%, 99% CL limits. 
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FIG. 6. Allowed (tan^ 623, A?tt,32) regions for different tan^ ^13 values, for the FINKS combina- 
tion of atmospheric neutrino data. The regions refer to 90, 95 and 99% CL. The best-fit point in 
the three parameter space is denoted as a star. 
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FIG. 7. Allowed three-neutrino {tan^ 623, /S.m'^2) regions for different tan^ ^13 values, for the 
CONT-UNBIN combination of atmospheric neutrino data. The regions refer to 90, 95 and 99% 
CL. The best-fit point is denoted as a star. 
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FIG. 8. Allowed three-neutrino (tan^ 023; A'7^32) regions for different tan^ 6*13 values, for the 
CONT-BIN combination of atmospheric neutrino events The regions refer to 90, 95 and 99% CL. 
The best-fit point is denoted as a star. 
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FIG. 9. Zenith-angle distributions for the Super-Kamiokande e-hke (left panels) and ^u-like 
(right-panels) events, both in the sub-GeV (upper panels) and niulti-GeV (lower panels) energy 
range. The thick solid line is the expected distribution in the SM. The thin full line is the pre- 
diction for the overah best-fit point of ALL-ATM data tan^ 613 = 0.025, Am§2 = 3.3 x lO'^ 
eV^ and tan^ ^23 = 1-6. The dashed (dotted) histogram correspond to the distributions for 
Am|2 = 3.3 (2.85) x lO^^ eV^, tan^ 6*23 = 3.0 (3.1) and tan^ 6*13 = 0.33 (0.54) which are allowed at 
90 (99)% CL. 
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FIG. 10. 90, 95 and 99% three-neutrino allowed regions in the (tan^ 023, ^77132) plane for 
different tan^ ^13 values, for the combination of UP-^ events induced by atmospheric neutrinos 
from Super-Kamiokande (through-going and stopping) and MACRO (through-going only). The 
best-fit point is denoted as a star. 
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FIG. 11. Zenith-angle distributions for upward-going niuon events in Super-Kamiokande and 
MACRO. The thick sohd Une is the expected distribution in the SM. The thin full line is the pre- 
diction for the overah best-fit point of ALL-ATM data tan^ 613 = 0.025, Am§2 = 3.3 x 10"^ 
eV^ and tan^ ^23 = 1-6. The dashed (dotted) histogram correspond to the distributions for 
Am§2 = 3.3 (2.85) x lO'^ eV^, tan^ 6*23 = 3.0 (3.1) and tan^ 613 = 0.33 (0.54) which are allowed at 
90 (99)% CL. 
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FIG. 12. Allowed (tan^ 023, ^77132) regions for different tan^ 6*13 values, for the combination of 
SK atmospheric neutrino events. The regions refer to 90, 95 and 99% CL. The best-fit point is 
denoted as a star. 



49 



10 



-2- 



-z- 



10 



-4 
,-^10 

(D _2|- 
^^10 

E -^ 
< 10 ^ 



10 



-4 



10 



-2- 



10 



-3- 



10 



-A 



::j.-_S. 



tan=iJt3=0. ::tan^iJij=0,026:: iarf'$^j=QAO: 




I I llllll| 1 I I llll j= I I I lllll| 1 I I llll l z 1 I I lllll| 1 I I I I I 1 ^ 

tan''i?i3=0.17!! tan^i?^a=0.26:! iov^%3=Q.5^ ! 




tan=i?ij=0.39 




J I 





I I M ill — I I iiiii jz I I I//III — I I iiiii iz / I m ill — I I iiiii ^ 



J I 





tan^tfij=0.45 f iarf4j3=Q35 : 



J I 



10 ^ 1 



10 1 10 



1 10 



FIG. 13. 90, 95 and 99% CL three-neutrino allowed regions in (tan^^23, ^1^32) ^^^ different 
tan^ ^13 values, for the combination of ALL- ATM neutrino data. The best-fit point is denoted as 
a star. 
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The dotted horizontal lines correspond to the 90%, 99% CL limits for 3 d.o.f. 
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FIG. 16. Allowed (tan^023) tan^^ia) regions for different A.m'^2 values, for the combination 
ALL- ATM of atmospheric neutrino data. 
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FIG. 17. Same as Fig. 8, but including also the CHOOZ result. The best-fit point is denoted 
as a star. 
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FIG. 18. Same as Fig. 13, but including also the CHOOZ result. The best-fit point is denoted 
as a star. 
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FIG. 19. Same as Fig. 16, but including also the CHOOZ result. 
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FIG. 21. Allowed regions in ATTi-g^ and tan^ ^12 from the global analysis of solar, atmospheric 
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